Ehrenfest relations at the glass transition: solution to an old paradox 
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In order to find out whether there exists a thermodynamic 
description of the glass phase, the Ehrenfest relations along 
the glass transition line are reconsidered. It is explained that 
the one involving the compressibility is always satisfied, and 
that the one involving the specific heat is principally incorrect. 

Thermodynamical relations are presented for non-ergodic 
systems with a one- level tree in phase space. They are derived 
for a spin glass model, checked for other models, and expected 
to apply, e.g., to glass forming liquids. The second Ehrenfest 
relation gets a contribution from the configurational entropy. 
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The glass transition is a dynamical freezing transition, 
that occurs when a liquid is supercooled. The transition 
is smeared, but becomes the sharper the slower one cools. 
For ideal, adiabatic cooling the transition will be sharp 
and occurs at the Kauzmann temperature Tk- 

Experimentally it is known that second derivatives of 
the free energy, the specific heat, the compressibility and 
the thermal expansivity, make a (smeared) jump from 
their liquid values to smaller values in the glass. Since 
many decades in time are involved, one might therefore 
wonder whether it can be described as a (smeared) sec- 
ond order phase transition. This idea has been put for- 
ward by Gibbs, DiMarzio and Adam. As thermodynam- 
ics amounts to system-independent laws, the approach 
leads to sine-qua-non relations along the glass transition 
line p{T). They are the Ehrenfest relations 
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where AA = Auquid - Agiass for any A. From exper- 
iments it was concluded that the first relation is usu- 
ally violated, while the second is closely satisfied in most 
cases, but not in all. (For reviews see |^). The 
Prigogine-Defay ratio 
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should be equal to 11 = 1. Experimental values are typ- 
ically found, however, in the range 2 < IT < 5. It is 
generally believed that 11 = 1 is a strict lower bound. 

In an extension of the theory one assumes that at the 
transition a number of (unspecified) internal variables Zi 



freeze in, and that the configurational entropy is constant 
along the transition line. This modifies (|l|) and (|2|) but 
keeps n = 1. Q These negative results have prevented 
further development of a thermodynamic approach. 

In this work we first explain that the first Ehrenfest 
relation is automatically satisfied (already in the reported 
measurements). By the same reasoning we shall conclude 
that the second relation must be incorrect. Analyzing 
model systems we shall then derive an extra contribution 
that arises from the configurational entropy. 

Before discussing the meaning of the Ehrenfest rela- 
tions, we first have to define the experiment, or, better 
said, the set of experiments, to be performed. Let us 
consider for definiteness cooling of a glass forming liq- 
uid at a fixed pressure pi and cooling rate Q = —dT /At. 
Starting from a high temperature one measures the spe- 
cific volume V{T]pi). It is linear at large T and at low 
T and has a smooth crossover between these behaviors. 
This happens near the freezing temperature Ti. Let us 
then repeat the cooling experiment at a large set of dif- 
ferent pressures p^, with moderate steps Pi — Pi-i. This 
will lead to a set of freezing temperatures T^, which de- 
fine a smooth freezing line p{T). The location of this 
line is by no means universal. It is defined by our set 
of measurements, here the set of pi's and their common 
value Q of the cooling rate. Different smooth sets of ex- 
periments may involve: a different cooling rate Q'; non- 
uniform cooling by letting Q — > Qi depend smoothly on 
Pi] non- linear cooling; or cooling where also p changes in 
time. All these sets of experiments will in principle lead 
to different transition curves p{T). In practice this means 
that old works in literature, where the cooling procedure 
has not been specified, are not reproducible. Likewise, 
computer experiments, with their extremely high cooling 
rates, cannot yield realistic glass transition temperatures. 

To test the first Ehrenfest relation (|^) one needs k, = 
~dh\V/dp\T, which is difficult to determine from cool- 
ing curves at two consecutive pressures. Therefore it has 
become standard to measure Ki in the glass phase by 
cooling at pi down below Ti and then making small pres- 
sure variations. |^ Such procedures, however, lead to a 
determination of dp /AT from an experiment at (essen- 
tially) pi only, a contradictio inter terminis! No exper- 
iment at one pressure can fix the slope of the transition 
line, because that depends on the conditions under which 
the set of experiments will be performed. ||^ The hope 
that the compressibility could be obtained by small pres- 
sure variations is frustrated by the history dependence 
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of the glassy state. A closely related phenomenon is 
known from experiments on spin glasses: the short-time 
("zero-field cooled") susceptibihty xzfc = (1 — qEA)/T 
is lower than the long-time ("field cooled") susceptibil- 
ity x_Fc=(l — /o dxg(x))/T, where qea is the Edwards- 
Anderson order parameter and where q{x) < Qea is the 
Parisi order parameter function. In the glass the short- 
time value of K, (measured by small pressure steps |^) 
will also be too low, yielding the observed too large Ak. 

The correct procedure is obvious and comes from the 
meaning of the Ehrenfest relation. The continuity of the 
specific volume can be considered at two glass transi- 
tion points (Ti,pi) and {T2,P2), where AV = 0. One 
may thus write AV{Ti,pi) - AV{Ti,p2) = AV{T2,P2)- 
AV{Ti,p2). (The terms at {Ti,p2) do not vanish). Di- 
viding by pi — p2 and taking the limit p2 —>■ Pi this, of 
course, leads to eq. (|l|). However, we must use on both 
sides of the equality the same values of pi —p2] it should 
not be taken "infinitesimal" on the left hand side and 
only "small" on the right hand side. In practice the only 
way is to determine k from the curves V{T;pi). Modern 
computer graphics allows to fit all the experimental data 
above and below the transition regions to high- and low- 
T surfaces in V — p — T space. Using all data should 
lead to reasonable fits. The intersection line of the sur- 
faces will satisfy eq. (|^) . This approach thus explains the 
old paradox: If properly interpreted, the first Ehrenfest 
relation is satisfied automatically! As there is no sec- 
ond procedure to produce the same glassy state, there is 
hardly a point in testing (|l]) experimentally. 

By the same token the second Ehrenfest relation 
(H), relating dp/dT to measurements at one pressure 
only, cannot be correct! Likewise, the Maxwell relation 
dU/dp + pdV/dp — —TdV/dT must be violated in the 
glass. We now show how unexpected behavior of the 
configurational entropy modifies it. Below we discuss the 
derivation for spin glasses, and then extend it to the glass 
transition in the hypernetted chain approximation. We 
believe that these relations are very general, and first 
formulate them for glass forming liquids. 

For glassy systems the entropy consists of two terms. 
Sint is the internal entropy, related to the glassy state 
the system condenses into; it lies well below the liquid 
entropy. Ic is the configurational entropy due to the 
number of equivalent glassy states. It is extensive at 
dynamical transitions, and becomes sub-extensive only 
for ideal adiabatic cooling. This part of the entropy 
is "lost" in the glass transition region. The quantity 
TdSint/dT\p will generally be smaller than the specific 
heat Cp = d{U + pV)/dT\p. In spin glass models the 
configurational entropy contributes to the free energy as 
— (T/x)Xc, where x is the weight (0 < a; < 1) occurring 
as breakpoint of the one-step replica symmetry breaking 
Parisi order parameter function, x also shows up dynam- 
ically as the factor by which the fluctuation-dissipation 



theorem is broken at long times. = T/x can be 
considered as effective temperature at which processes 
related to the configurational entropy thermalize. ^ 
It was shown that Cp = TdS^nt/dT + TMc/dT. 
We are now in the position to take the first derivatives 
of the free enthalpy G ^ U + pV - TSmt - T^Ic- This 
yields 



^-V-^I (4) 

dp~^ dp ' ^ ' 



Along the transition line G, Suq = Sint + 1c and V are 
continuous wrt to the liquid. The standard assumption 
that the first derivatives of G are also continuous, is seen 
to be incorrect |^: the terms involving arc nontriv- 
ial in the glass {dT^/dT < 1, dTe/dp ^ 0). The finite 
diff'erence in slopes discussed here is due to unexpected 
behavior of the configurational entropy, [|j neglected so 
far. It leads to the modified Maxwell relation 

\dV pdV dV _ dT^ die dT^dl^ 

Along the freezing line one has Te{T,p{T)) = T and 
one may define the total derivative d/dT — d/dT + 
{dp/dT)d/dp. Eq. (^ does not violate the balance: 
dAG/dT = since dTjdT = 1. Let us now con- 
sider eq. (H) and subtract the values on the liquid- 
side. Multiplying by dp/dT and using dTg/dT = 1, 
dAJ7/dT=dAy/dT= we obtain from (|) the modifi- 
cation of the second Ehrenfest relation (c.f. eq. 



ACp 
TV 



Aa 



dp 
dT 



(1 



STe. die 



dT ' VdT 



(6) 



This relation indeed connects dp/dT with another deriva- 
tive along the transition line, namely that of the configu- 
rational entropy. This term originates from the difference 
in slopes of the liquid and glass free enthalpies. The fac- 
tor 1 — dTe/dT > is a nontrivial weight. Since the first 
Ehrenfest relation is satisfied, measurement of n is not 
needed for the Prigogine-Defay ratio: 11 = IldT/dp with 
n = ACp/ [TVAa) determined at one pressure. 11 will 
be less than unity when dl^/dT < 0. 

We now consider the data of Rehage and Oels for the 
glass transition of atactic polystyrene. ||l^ For cooling 
at a speed of 18 K/ hour at p = 1 bar they report: T — 
361 K, ACp/V = 0.30 J/gK, Aa = 3.510-4 cm^/gK, 
dp/dT = 0.31 bar/K, Ak = 1.6 10"^ cm^/gbar. This 
was reported to yield 11 — 1.06 « 1.0, and a violation of 
the first Ehrenfest relation. This violation has, however, 
already been traced back to the way how k was measured. 
Our Prigogine-Defay ratio 11 = IldT/dp = 0.77 is less 
than unity. The last term in eq. (^ is negative and brings 
23% of the value for the slope dp/dT, a large effect. 

We have discussed how the configurational entropy 
modifies an Ehrenfest relation. This effect should be 
stronger for first order glass transitions, which occur, for 
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instance, in water, In the p-spin model (see eq. (||)) 
this happens when the transversal field T exceeds a criti- 
cal value. Using that along the dynamical transition 
line AG = dAG/dT and eq. (|) we obtain 

AV^ = ^{AU+ pAV) + (| - (7) 

which deviates by the Xc term from the static Clausius- 
Clapeyron equation. Note that x = T/T^ is below unity. 

Let us now give the theoretical background of our 
relations. They have initially been derived within a 
spherical p-sp'm interaction spin glass. For a system 
of m-component spherical spins Sf {i = l,---,N, c — 
l,---,m), satisfying J2ic^i^ ~ consider the 

Hamiltonian in a transversal field 
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The independent Gaussian quenched random couplings 
have average zero and variance J^p!/2A^''~^. The system 
has a multitude of states a — 1, • • • , A/", each with its 
own free energy Fa that is a local minimum of a known 
TAP-free energy functional. The replica calculation with 
one-step replica symmetry breaking involves parameters 
/i, qd^ qi, and x, leading to the free energy pTI] 
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where = 1 — x. ^, qd, and qi are determined by op- 
timizing F. It was recently pointed out by us that the 
value of X is related to the time scale at which the system 
is considered, jl^ Indeed, setting dF/ dx = (leading to 
eq. (^o|) with r] ^ rjst < ^ independent of T and F) yields 
the static phase transition at the Kauzmann temperature 
Tk, related to the longest time scale. On the other hand, 
the marginality condition (eq. ( |l0|) with ?; = 1) describes 
algebraic time scales, at which a transition occurs at a 
higher temperature Tj^. This is reproduced by the mode- 
coupling equations, and is comparable to the sharp crit- 
ical temperature occurring in mode-coupling equations 
for glasses. As this transition is absent in practice, we 
have considered the system at exponential time scales 
t = ioexp(iVr). At given r barriers with free energy 
height less than NTt can be surpassed, r/ parametrizes 
the value of the free energies of the states a. As time 
evolves, the dominant lowest reached free energy Fmin(t) 
has parameter ri{t). When, at fixed field Fi, also the 
temperature T{t) is slowly lowered, we can eliminate t to 
obtain a function ?y(T; Fi). Cooling trajectories at a large 
set of fields F, define a set of experiments. If the set is 
"smooth" it will lead to a smooth function 77(T;F). This 
function should follow from solving the dynamical equa- 
tions. We shall not do that, but remain on a quasi-static 



level, where the information of the cooling dynamics is 
coded in the function ri(T]T). 

In our present analysis a freezing transition occurs 
when the temperature, below which the dominant lowest 
reached state will freeze, is equal to the actual tempera- 
ture. We assume that we can still describe the situation 
by the Gibbs weight, which is the case when no relevant 
parts of phase have become inaccessible. 

The free energy of the TAP states can be characterized 
by a parameter rj (rjst < ?7 < 1), that enters the condition 
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In solving the saddle point equations for /i, g^, and gi, 
the above relation fixes x = (j> — 1 — rj) {qd ~ qi) /vQi ■ One 
can calculate all quantities of interest. We have verified 
the following relations for i^(r, re(r, F), F): 
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This implies the modified Maxwell relation 
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-F, U , Sint, Ic, M and Te only depend on the value of 
rj in the point (T, F) . Their temperature derivatives also 
depend on drj/dT, while their field derivatives involve 
drj/dp, the measure of variation between experiments at 
different fields. 

The second Ehrcnfest relation can be rederived by mul- 
tiplying eq. (ph by (l/7V)dF/dT. This yields generally 
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while Aa — A^dF/dT is again satisfied, where a = 
-{dM/dT)/N and x = {dM/dT)/N. In comparing with 
(|) one should keep in mind that dT^/dT = 1. Eq (|l|) 
also applies to the unconventional case where one starts 
cooling adiabatically below Tk- At freezing one then has 
le — (Me/ AT = 0, but the last term is non-zero. 

dlc/dT does not depend on drj/dT but only on dT/dT. 
It holds that AG, Aa and Ax are proportional to 1 — 
dTe/dT (since near the transition Usg — UpM ^ x — \). 

The relations (|l^^5|) are expected to be universal for 
dynamical glassy transitions with extensive configura- 
tional entropy. We have considered three other cases: 

1) In the case of a longitudinal field (F^S'f — > 
Hj2^i) the same equations are satisfied, see also JT2[ . 

2) The hypernetted chain equation of fiuids is an 
approximate non-linear integral equation for the pair- 
correlation function, Mezard and Parisi M pointed 
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out that in a certain region it has many solutions, describ- 
ing a glass phase. By weakly coupling different copies 
(replica's) of the system, they introduced a replica cal- 
culus. The main difference with the above spin glass is 
that the spin-spin overlap is replaced by the pair corre- 
lation function. The static transition again follows from 
the relation dF/dx = 0. Dynamically this relation is 
not satisfied, and there is an extensive configurational 
entropy. We now can make the same assumptions as in 
the above spin glass and will rederive eq. (^). 

3) Recently we have introduced a model of a directed 
polymer on a square lattice with a correlated random 
potential, consisting of randomly located parallel ridges 
(repulsive potentials). The polymer prefers to lie in 
broad lanes (width £) in between ridges, of which there 
occur a lot when the transversal width scales as W = 
exp(AL"'^/^) in the parallel width L. The free energy reads 
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where /s(T) is an uninteresting bulk free energy den- 
sity, r is the interface stiffness, exp(— /i) is the chance for 
having no ridge at a given height, and j/ is a chemical 
potential favoring (y > 0) or disfavoring < 0) wide 
lanes, and Xc = log W — ^£ is the complexity. At some 
Tk, where j{T;u) = {^l/X)[TT{T)/iT^l - i^)]^^^ equals 
7(T;j/) = 1, there is occurs a static "Kauzmann" transi- 
tion from a high-temperature phase, where the polymer 
lies in the broadest lane {£ = imax = AL-^/^/^), to a low 
temperature phase, where it spends most of time in a set 
of narrower lanes {£* = jimax with 7 < 1). 

Starting from a large set of uniformly distributed in- 
dependent polymers, we are interested in the dynamical 
(short time) regime A, where (M) is valid with £ increas- 
ing logarithmically with time. ||l5tl In order to make con- 
tact with previous theory, we replace TXc in (|6|) by TeTc 
where Te{t)=LTT{T)/{fie{t)^) + ly/fi. A minimum will 
then occur a,t £ — £{t). In the polymer model one has a 
reversed role of heating and cooling. jTsj A set of experi- 
ments can be introduced by specifying smoothly related 
heating trajectories at a large number of ViS. For each 
of them a dynamical phase transition can occur at any 
T < Txii^i), where the dominant width £(t) reached so 
far, equals its freezing value £*{T{t)). This transition is 
very similar to the above ones. Eqs. ( pT[]l5| ) are satisfied 
with r ^ I/. As before, Tc is large at the transition. The 
second term in the r.h.s. of eq. ( |l5| ) is finite, while the 
last one vanishes. A related transition can occur at any 
T > Tk. Then the transition line is £{t) — £max, where 
Xc = 0. Now the last term in ( [Tsl ) is non- vanishing. 

In conclusion, we have pointed out that the present 
understanding of the Ehrcnfest relations is incorrect. We 
have explained that the first one is satisfied automatically 
and that the second one must be modified. From quasi- 
static model calculations we have shown that it gets an 
extra contribution from the configurational entropy. This 



explains that the Prigogine-Defay ratio is smaller than 
unity in the experiments or ref. [Tc| ] 

We have also presented the generalization of the stan- 
dard thermodynamical laws to non-ergodic situations 
with a one-level tree in phase space. This is given in 
eqs. (pTl)-(p^, where Tg is an effective temperature, that 
depends slowly on time. Along the transition line the 
modified Maxwell relation eq. ( p^ leads to a new form 
(p|,p^ of the second Ehrenfest relation. We have checked 
the predictions in several model systems. 

Note added: we have checked our Ehrenfest relations 
for cooling in the backgammon model, which has no dis- 
order but entropic barriers. |17[] 
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